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Abstract 

In this paper, we study the steady solutions of Euler-Poisson equations in 
bounded domains with prescribed angular velocity. This models a rotating New- 
tonian star consisting of a compressible perfect fluid with given equation of state 
P = p"' . When the domain is a ball and the angular velocity is constant, we 
obtain both existence and non-existence theorems, depending on the adiabatic 
gas constant 7. In addition we obtain some interesting properties of the solu- 
tions; e.g., monotonicity of the radius of the star with both angular velocity and 
central density. We also prove that the radius of a rotating spherically symmet- 
ric star, with given constant angular velocity and constant entropy, is uniformly 
bounded independent of the central density . This is physically striking and in 
sharp contrast to the case of the nonrotating star. For general domains and 
variable angular velocities, both an existence result for the isentropic equations 
of state and non-existence result for the non-isentropic equation of state are also 
obtained. 

1 Introduction 

The purpose of this paper is to understand the structure of a rotating star. In the 
theory of General Relativity, for a non-rotating star consisting of a perfect fluid, the 
exterior and interior solutions are very well understood. The exterior solution of the 
star, that is, the gravitational field (space time metric) outside the star, is given by 
the well-known Schwarzschild solution, and the interior solution, which matches the 
exterior Schwarzschild solution on the boundary of the star, is obtained by solving the 
*This paper has appeared in Arch. Rational Mech. Anal. 173 (2004) 345-377. 
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Tolman-Oppenheimer-Volkoff equations ^8]. For a rotating star, the exterior solution 
of the gravitational field is the celebrated Kerr solution jl] . However, no corresponding 
interior solution of Einstein's equation is known which matches the Kerr solution at the 
boundary of the rotating star. Actually, even in the Newtonian case, unlike the case 
of the non-rotating star which has been extensively studied by Ritter, Lane, Emden, 
Kelvin and Chandrasekhar; (see j3] ; and also |H] and ^3] for niore recent results in the 
case of non- rotating star), there have been very few results on the rotating star if the 
star consists of a compressible perfect fiuid. This paper is devoted to the investigation 
of the steady solutions of a compressible perfect fiuid rotating star in the Newtonian 
case with prescribed angular velocity, in a bounded domain, and with zero density on 
the boundary chosen to match the exterior of the star. 

The motion of the compressible perfect fiuid with self-gravitation is modelled by the 
following Euler-Poisson equations; cf [H], 

pt + div^(pv) = 0, 
pwt + (pv ■ V)v + VP + pV$ = 0, 
5i + V ■ V5 = 0, 
A<l> = A-nGp, 

where p, v, P, S and $ denote the density, velocity, pressure, entropy and gravitational 
potential, respectively. Here (t, x) G x denotes the time and space variables and 
X = {xi,X2,X3). For simplicity, we assume the pressure satisfies the following equation 
of state 

P = P(p, S) = eV, (1-2) 

where 7 > 1 is the adiabatic exponent. 

Suppose the star rotates about the xs-axis; we are interested in finding an axi- 
symmetric solution (p, v, 5*, t) = (p, v, S, ^){ri{x),X3,t) of (jl.lj) with pre- 

scribed time-independent angular velocity ^{rj), where 



T] = 'q(x) = \J x\ + 



In this case the velocity field is given by v = {—X2^{rj) , xiVLirf) , 0), and thus div^(pv) = 
and V ■ VS = 0. By (l.l)i and (1.1)3, = and St = 0. Thus, the solution 
(p, V, 5*, $) of (1.1) is time- independent and satisfies the following system of equations 

VP(p, ^) + pV$ - pV J(r^) = 0, ^^^^ 
A$ = AnGp, 
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where 

rv 



J{r]) = / sn'{s)ds, (1.4) 







and the entropy S" is a given bounded function. 

In for an isentropic fluid, i.e., the entropy S = constant (without loss of general- 
ity, S is assumed to be zero in this case), Auchmuty and Beals considered system ()1.3p 
in all of R^, with prescribed total mass M, and gravitational potential $ given by 



$(a;) = -G [ T^^dy =: -Bp{x) 



The problem then reduces to finding the solution of the following equation, 

V(^-i?p- J(r/))(x) = 0, (1.5) 
7 - 1 

where p > 0. They formulated this as a variational problem; namely, minimize 

E{p)= / {-^--p-Bp-pJ{r,)){x)dx, (1.6) 

in the class 

VTm = {p > : / p{x)dx = M}. (1.7) 
By assuming that the angular velocity Q satisfies the following decay properties 

J(+oo) < +oo,J e C^[0, + 00 , 

?7(J(+oo) - J{r])) ^0, as 1] +oo, (1.8) 

and that the adiabatic exponent 7 satisfies 

7 > 4/3, (1.9) 

Auchmuty and Beals proved the existence of a minimizer of the functional E{p) in the 
class of functions (|1.7|) . Moreover, this minimizer has compact support and satisfies 
equation (|1.5p wherever it is positive. The shape of the free boundary which separates 
the vacuum and fluid was investigated in |3| for the Auchmuty-Beals solutions. The case 
of an isentropic uniformly rotating star, (i.e., the angular velocity Q is constant), was 
discussed by Li in (^I]); he proved the existence of a minimizer of the functional (1.6) 
in the class Wm , under the assumption 7 > 4/3. The diameter of the support of the 
density p was studied in [H] for the solution obtained in [TTj. In the proof of the above 
results, the prescribed total mass serves as a constraint on these variational problems. 
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and without this constraint, it is not clear that the minimizer of E{p) exists. In 
the angular velocity is prescribed in the entire space R'^ (even in the vacuum region), 
and is assumed to satisfy the decay property ()1.8j) . In [2j and [TT], 7 is required to be 
greater than 4/3, so a natural question one can ask is what happens when 1 < 7 < 4/3 
? Another issue is that all of the above mentioned results are for isentropic fluids, so 
one can also ask what happens if the fluid is non-isentropic. The purpose of this paper 
is to address these, and other issues. 



2 Statement of Results. 

We are interested in the solution of ()1.3|1 in a bounded domain D in R^, where p{x) > 
for x G -D and p{x) = if x G dD. From the first equation in system (1.3), we have 

lvP = V(J-$) (2.1) 
P 

for X G -D. Thus, by the second equation of (1.3), we obtain 

divi-VP) = A(J -<^) = AJ -inGp (2.2) 
P 



for X E D. Set 

7 



w 



7-1 

Then it is easy to verify, using (1.2), that 



IvP = e^'^S/w. (2.4) 
P 

Substituting this in ()2.2|1 . we obtain the following elliptic equation 

div^e'^^Vw) + Ke-'^^w" - 2n{r]){n{r]) + r]n'{r])) = 0, (2.5) 

where 

q= -, a = -, (2.6) 

7-1 7 

K = AtxG (^^^j ^ ^ ) ci^d for simplicity, we can normalize K to make K = 1. We seek 
solutions of (|2.5|) satisfying 

p{x) > 0, xeD, p{x) = 0, a; G dD, (2.7) 

or equivalently 

w{x) > 0, xeD, w{x) = 0, X G dD. (2.8) 
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In this paper, we only consider the case when 

1<7<2, 

because when 7>2,0<g<l, and equation ()2.5|1 becomes subhnear, and this 
situation was studied completely in Moreover, equation ()2.5j) is linear if 7 = 2, 
and there is a complete theory for linear elliptic equations (cf. jlHl)- 

We first consider the case when the angular velocity Q{ri) = Q = constant, the 
entropy S{x) is spherically symmetric ( i.e., S{x) = S{r), r = \x\ = \/ a;|), 
and the domain D is a ball i?ij(0). In this case, we look for the spherically symmetric 
solutions of problem ()2.5|) and ()2.8|) . i.e. , w[x) = w{r), r = \x\ = >y xf + xl + x|. 
Then w{r) satisfies the following equation 

w"(r) + (- + aS'(r))w'(r) + e~"^('^)(e-"^(^')w« - fx) = 0, (2.9) 
r 

and boundary conditions 

w'(0) = 0, w{R) = 0, (2.10) 

where we have set 

cr = 2n'^. (2.11) 

We want to use the "time-map" method used by Smoller and Wasserman (cf JZI). For 
this purpose, we consider the problem ()2.9|1 with initial data 

w{0) =p>0, w'{0) = 0, (2.12) 

where p is a free parameter, and let w{r,p,a) be the solution of this problem. We 
define the " time-map" p — > R{p, a) , by 

R{p,(t) =m{{R\ R> 0, w{r,p,a) > if < r < R, and w{R,p,a) = 0}, (2.13) 

so R{p, 0") is the first "time" at which ti; is (we will write R{p, a) = 00 if w{r, p, a) > 
for all r > 0. ) Thus R{p, a) is the radius of a rotating star with given central density 

(cf. dOD) 

p(0) = (I^p)^e-^W/^. 

7 

and angular velocity 

n = 

In order to state our first theorem, we introduce the following notation. We set 

^ = sup5(r),^=infS'(r). (2.14) 

r>0 '">0 
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We assume 

-oo<S<S<oo. (2.15) 

We define the constant b by 

(1 + 2e"^)[10(l + g)]5(5 - 
Obviously 6>0ifl<g<5. Our first result is the following theorem. 

Theorem 2.1. Assume that D is a ball in and that Q is a nonzero constant. 

1) If 1 < 7 < 6/5 an(^ S'{r) >Oforr>0, then 

w{r,p, a) > 0, 

for every r > and p > . 

2) If 4/3 < 7 < 2, then 

i?(p,a)<+oo, (2.17) 
for sufficiently large p, provided the entropy satisfies the following conditions 

aS'{r)e''^z' > -^z, (2.18) 



for < r < Tc/y/b, wh 



ere 



, , sin(v6r) 
Ar) = (2.19) 
vbr 



and supo<r<i4.^ I'S"(?")| is sufficiently small. 

— Vb 

3) If 6/5 < 7 < 2, assume that the conditions in 2) hold, and assume too that the 
entropy satisfies the following condition 

S'{r) < 0, for < r < (2.20) 

y/b 

where b is given by \2.1t^) . Then 

R{p, a) <+oo, (2.21) 

for sufficiently large p. 

4) For the solutions w{r, p, a) =: w{r) in 2) and 3), we have the following estimates 
on the mass M{r) = 4:7iT'^p{T)dT, and the average density p{r) = yp^^ (o) "^(^) 
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ball Br{0), (here p = 



1 1/(7-1) 



IS the density (cf. (2.3)) and Vol 5,,(0) = Airr^/^): 



M{r) > 



2r^ 
3G 



(2.22) 



p(r) > 



(2.23) 



2nG' 



forr < R{p, a), where G is the Newtonian gravitational constant. 

Remark 1. The conditions imposed on the entropy S{r) in Theorem 2.1 are automati- 
cally satisfied in the isentropic case, i.e., S{r) = constant. 

Remark 2. The existence of positive solutions for the equation of isentropic non-rotating 
star is standard (cf. [5], [12] and [15]). In [8], for the case of a non-rotating star {Q = 0), 
the authors first consider non-isentropic equations of state, and an existence theorem 
is proved for a bounded domain D, under the assumption that the entropy function S 
satisfies, for some constant a > 0, 



for all test functions G Hq{D). Our condition (2.18) is an explicit condition on S 
and does not involve the second derivative of S, for the case when domain Z) is a ball. 

The proof this theorem, as well as the other theorems stated in this section will be 
given in sections 3-5. 

In the case of isentropic fiuids (i.e., S = constant), we can obtain further results on 
the qualitative properties of the solutions if the angular velocity fl{ri) = Q = costant 
and the domain D is a ball. Without loss of generality, we may assume 



In this case, it follows from the celebrated Gidas, Ni and Nirenberg result (f^) that 
positive solutions to (|2.5|) and (j2.7p must be spherically symmetric. Substituting S = 
in ()2.9p . we obtain the following equation: 




S = 0. 



(2.24) 



w"(r) H — w'(r) -\- w'' — a = 0. 
r 

We consider the problem ()2.25|) with initial data 



(2.25) 



w{0) =p>0, w'{0) = 0. 



(2.26) 
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We will again use w{r,p, a) to denote the solution of the above problem. The following 
theorem gives some physically interesting properties of R{p,a), the radius of the star 
with central density (cf. (2.3) with 5 = 0) 

P(0) = (^^P)^. 

7 

Theorem 2.2. Assume that D is a ball in , Q is a nonzero constant and the entropy 
S is constant (we set S = for convenience). Then the following statements hold. 
1) If 6/5 < 7 < 2, there exists a constant po > depending only on 7 and a such that 



'5a(g+ 1)" 


1/9 


'l07^]2" 


7-1 


5 — g 




57 — 6 





R{p, a) < +00, if and only ii p> po. 
Moreover, we can estimate Pq from below, 

Po > 

2) If 6/5 < 7 < 2 , 
provided ai > a2 > 0, and 
ifP2 > Pi> Po- 



R{p, CTi) > R{p, 0-2) 
R{pi,(^) > R{p2,(y) 



(2.27) 

(2.28) 

(2.29) 
(2.30) 



As a corollary of this result, we have the following theorem, in which pq is the positive 
constant given in Theorem 12.21 

Theorem 2.3. Assume that D is a ball in R^, Q is a nonzero constant and the entropy 
S is constant (set 5* = for convenience) . If 6/5 < 7 < 2, then there exists a positive 
constant C depending only on po and Q such that 



7-2 



7-2 



R{p, a) < Cp^^ < Cp^'-''-'^ , 
for p > Pq. In particular, this implies 

R{p, cr) ^ 0, as p ^ 00. 
We now discuss the physical meaning of the above two theorems. 

Physical meaning of Theorem\ 



(2.31) 



(2.32) 
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Part 1) means, for 6/5 < 7 < 2 and cr = 2Q^ > , there exits a critical central 
density, determined by po, depending only on 7 and the angular velocity fl, such that 
the radius of star is finite if the central density is greater than or equal to this critical 
central density; otherwise the radius of the star is infinite. This is in sharp contrast 
to a non-rotating star, for which in this range of 7 the radius is finite no matter how 
small the central density is (cf. The reason for this is that for the rotating star, 
the central density must be large enough to provide sufficient gravitational attraction 
to balance the centrifugal force due to rotation. The physical meaning of part 2) is 
that, for a compressible fluid the radius of the star increases with increasing angular 
velocity for fixed central density, while the radius of star decreases with increasing cen- 
tral density for the fixed angular velocity. 

Physical meaning of Theorem \2JA 

If 6/5 < 7 < 2, in order to have the radius of the rotating star R{p, cr) to be finite, p 
has to be greater than or equal to Pq, and the constant po is completely determined by 
7 and the angular velocity f2, as stated in Theorem 2.1. Thus, ()2.3H1 in Theorem 12.31 
shows that there exists a (finite) upper bound on the radius of a rotating star, which 
is determined only by 7 and f2, and is independent of central density, for the radius 
of the rotating star to be finite. This is physically striking and completely different 
from non-rotating star, for which the radius can be arbitrarily large. Actually, for the 
non- rotating star, the central density can also be arbitrarily small (cf. [5J. Moreover, 
as we will show later, the radius of a non-rotating star is proportional to (p(0))^~, 
where p(0) is the central density. Since 7 < 2, the radius of a non-rotating star can 
be very large if the central density is very small. However, for the rotating star, the 
central density cannot be arbitrarily small, as indicated in part 1) of Theorem 12. 2j cf. 
(2.27). Moreover, (j2.32|) shows that the radius of a rotating star tends to zero if the 
central density tends to infinity. 

We now turn to the case of general bounded domain D C and variable angular 
velocity ^{rj). Without loss of generality, we may assume that G -D. For the existence 
of positive solutions to the boundary value problem ()2.5j) and ()2.7j) . we assume that 
the angular velocity VL{ri) is as a function of r] and satisfies the following condition: 

Condition A: Let 5/^^(0) = {x\ \x\ < Ri} be the larggest ball contained in D which 
is centered at origin. For x E D — Brj^{0), assume 

nmv^{v)]'<o, (2.33) 
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and is not identically zero; here, as before, r] = rj{x) = ^/ xf + x^, and " prime " denotes 
A 

dri ' 

Remark 3. Notice that ()2.33|) is only required for x E D — B^^iQ). Therefore the 
angular velocity ^{rj) can be any function when D is a ball. 

The following theorem is an existence theorem for a general domain D and variable 
angular velocity. 

Theorem 2.4. Assume 6/5 < 7 < 2, the entropy S is constant (set S* = for con- 
venience), VL{rj[x)) G C^^^{D) n C{D) for some 5, < 5 < 1 with dD smooth. Then 
there exists a constant R2 > depending only on Q and 7 such that if D G B^.^lO) and 
condition A holds, there exists a positive solution to i2.^} and \2. 7| ). Moreover, for this 
solution , the central density satisfies : 

p(0)>(^)^/(^-^)--^-/3, (2.34) 
7 07 — 

where P = max^.^^ \2fl{rj){fl{r]) + riQ'{ri)\. 

Notice that as in Theorem 1.3, the size of the domain D is uniformly bounded, where 
the bound depends only on 7 and fl. 

Let w{x) be the positive solution given in Theorem 12.41 Notice that since 5* = 
(2.3) implies p = {^^w)~ . We can obtain the potential function by virtue of (1.10) 
and (1.13). Actually, when S = 0, (2.4) and (2.5) imply V(w + $ - J) = 0. So we 
set $(a;) = {C — w + J{ri)){x) for x E D, where C is any constant. Then (p, <l>)(x) is 
a positive solution to the system (1.3) in the domain D. System (1.3) is derived by 
assuming that p is axi-symmetric. In fact, we have the following self consistency result. 

Theorem 2.5. Suppose the conditions in Theorem I.4 hold. If the domain D is axi- 
symmetric about X3- axis, then the solution w{x) in Theorem \2.4\ is also axi-symmetric 
about x^-axis. 

Our final result is the following nonexistence theorem; here v{x) denotes the unit 
outer normal vector on the boundary, and rj = rj[x) = \fx\^rx\\ 

Theorem 2.6. Suppose dD G and the domain D is bounded and star-shaped, i.e. 
X ■ i^ix) > for X G dD, and the entropy S satisfies the condition 

X ■ VS{x) > 0, for X ED. (2.35) 
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If 1 ^ 6/5 (resp. 7 < 6/5) and 



5A{r]) 
2 



> {resp. > 0) 



(2.36) 



forx G D, where A{r]) = 2Q{ri){riQ{ri)y , then there is no positive solution to \2. 5\) and 



Remark 4. Condition ()2.36|) is trivially satisfied when is a constant. 

Remark 5. In [8], for the case of a non-rotating star [Vt = 0), the authors prove a 
non-existence theorem for 1 < 7 < 6/5, under the assumption that the entropy S{x) 
satisfies the following three conditions, for all x G -D: 



In Theorem 2.6, we only require x ■ VS* > for all x G -D. 

We now make some remarks about the above theorems and their proofs. In the 
theory of second order elliptic equations, the existence of a positive solution to the 
boundary value problem 



where D G -R" is a bounded open set, has been extensively studied either by the 
Min-Max method of the Mountain Pass Lemma of Rabinowitz jT3] or the Topological 
Degree Method (Leray-Schauder degree) by P. L. Lions fT^, both under the condition 
/(x, 0) > 0. If /(x, 0) < 0, the difficulty in solving the boundary value problem ()2.37p is 
that the Harnack inequality is not applicable, and thus the above variational methods 
(Mountain Pass Lemma), and the Topological Degree Method do not work in this case. 
The first existence result of a positive solution for the case /(x, 0) < was given in 
([T7]) by using the time- map method when D is a ball and /(x, w) = f{w) (i.e., / does 
not depend on x explicitly), this corresponds to the case when both the entropy S and 
the angular velocity f2 are constant in equation ()2.5j) . For this type of equation, further 
analysis was given in The existence result of a spherically symmetric solution to 
()2.5|) and ()2.8|) generalizes the above mentioned results to the non-isentropic case. 

If D is a general domain, ^\ contains some existence results for the case f{x,w) = 
f{w) and / is sublinear in w. To the best of our knowledge. Theorem 12.41 is the first 
result of existence of positive solutions for the case that /(x, w) is superlinear in w and 
the domain D is different from a ball. 



X ■ > 0, AS > and V(AS) ■ x > 0. 



Aw + f{x,w) = 0,x e D 



(2.37) 
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The rest of this paper is organized as follows. In Section 3, we prove Theorem 2.1. 
Theorems 2.2 and 2.3 are proved in Section 4. Theorems 2.4 - 2.6 are proved in Section 
5. Further discussions are given in Section 6. 



3 Proof of Theorem 2.1. 

In this section, we shall prove Theorem 2.1. We define 

f{r, w) = l ' - ' (3.1) 

[ -0-, if w <0, 

and 

F{r, w)= f{r, z)dz. (3.2) 
Jo 

As in Section 2, we use the notation w{r,p, a) to denote the solution of the problem 

w"(r) + (- + aS'(r))w'(r) + e~"^f(r, w) = 0, (3.3) 
r 

and 

w'{0) = 0,w{0) =p> 0. (3.4) 

Let R{p, a) be as in ()2.13|) : i.e., R{p, a) the first point at which w is 0. 

Motivated by the celebrated Pohozaev identity (^14 ), we define the function G{r) 

by 

Gir) = /'e'^M!^ + F(r, w)) + Ve"^wti;'(r) + | [\^e''^ SWdt, r > 0. (3.5) 



It is easy to verify that 







G'(r) = 3r\F(r, w) - l^fl^^) + r^gU, w), (3.6) 

6 

where 

M^..= (3.7) 

0, w <0; 



dr 



thus 

G{r) 



l\3t'[Fit,w) - l^^l^]+t^g{t, w)}dt, r > 0. (3.8) 



We now prove that 

// 1 < 7 < 6/5 (q>5) and S'{r) >0 forr>0, then 



R{p, a) = +oo (3.9) 
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for any p > and cr > 0. 

The proof is by contradiction. Suppose R{p, cr) < +00, and let R = R{p, a)). Thus 

w{R) = 0, w{r) > ioT < r < R. (3.10) 

By the definition of G{r) (3.5), we have G{R) = ^^'^"^'^^'K(^)^ + | r^e'''^ S'{r)w'^dr. 
Applying this to (3.8), we obtain, by virtue of ()3.H) and ()3.2|1 . 



i?3e"^(-R)(w;'(i?))2 



a 



2 2 



R 



+ n / r'e''^S'{r)w'^dr 







= ^'^3r2 {[^^ " ^)e-"V^ - ^aw;] - r3^^5'(r)e""^ti;^+'| cir. (3.11) 

Since — | < for g > 5 and S"(r) > 0, p.llj) gives a contradiction. Thus ()3.9p 
holds so Part 1) in Theorem 2.1 is proved. 

Next, we give the proof of parts 2) and 3) of Theorem 2.1. First, it follows from 
(rO|l that 

(r^e^V)' + r^f{r, w) = 0. (3.12) 
Thus from ()3.12|) and (2.4), we obtain 

w'{r) = -^^^f t^f{t,w{t))dt, r>0. (3.13) 
1" Jo 

For fixed p > 0, since w{0) = p, by ()3.1|1 . ()3.2|1 and ()3.13j] . we have w'{r) < for 
small r, r > 0, if p is sufficiently large. We define ri to be the point such that 

w{ri) = 2p/3, w{r) > 2p/3 for < r < ri. (3-14) 

Then we have the following lemma which estimates ri and G{ri) in terms of p. 

Lemma 3.1. For sufficiently large p , if 1 < q < 5, (6/5 < 7 < 2j, then there are 
positive constants C\, 02 and C3 independent of p such that 

c,p-(i-^y^ < n < C2p-('?-^)/^ (3.15) 

and 

G(ri) > C3p(^-«)/^ (3.16) 
prowded supo<r<ri l*S"(r)| is sufficiently small. 
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Proof. By the definition of ri (cf. ()3.14j) . we liave 

2p/3<w{r)<p, 0<r<ri. (3.17) 
Tlius, from dSIU, (j3T3|l and (jTTTjl . 

/•ri -aSir) rr 

2p/3 = w;(ri) =p- j —j^ / t2(e-°^(*)«;^(t) - a)dtdr 



>p- fie-"^p^ - a)dtdr 

Jo Jo 



2 

= P--(e"'"V-e-"^a), (3.18) 
6 

if p is sufficiently large, where S_ = info<r<+oo 5'(r), and by (2.15) S_ > — oo. From 
dSmi), we get 

2 

^(e-^°V - e-"^a) > p. (3.19) 

This implies 

for some positive constant ci, if p is sufficiently large. Similar to the argument in ()3.18|) . 
we can obtain, 

2p/3 < p - d{e-'"S{2p/3y - e-"^a), (3.20) 
o 

where S = supo<r<+oo 'S'(r), and S < oo. So from ()3.2()|1 . 

for some positive constant C2, if p is sufficiently large. This proves (3.15). To prove 
()3.16|) . we have, in view of ()3.8|) . 

Gin) = r 3r^ I [— - , " , rgVr)]e-°^w'?+^ - -au^l rfr. (3.21) 

^ ^ 7o Vq + l 6 3(g+l) ^ 6 J ^ ^ 

Since 1 < g < 5, (g + 1 < 6), we have 

1 1 q; 

^ - - ^^TTTTT^^V) > ^3, 0<r<ri, (3.22) 
g+1 6 3(g+l) 

if suPq<^<^^ I'S"(?")| is sufficiently small, where C3 is some positive constant independent 
of p. By virtue of (jlTTll . (HT^ and (jT^ . we have, 

G(ri) > (^C3e-^(«;(ri))'^+^ - ^ap 

= rl fc3e-"^(2p/3)''+^ - ^ap^ . (3.23) 



6 

Therefore, if p is sufficiently large, ()3.23|) implies ()3.16p . This proves Lemma 3.1. □ 
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By (jSUl), we have 

/(r, w) > 0, for w > [e^^a]^/". (3.24) 

Therefore, in view of ()3.13p . we have w'{r) < for small r if p is sufficiently large, so 
w{r) decreases for small r. Using ()3.13|) and ()3.24|) . we can see that w'{r) < and thus 
w{r) decreases as long as w > [e^'^cr]^/'^. Because of this, we can define T =: T(p) as 
the point such that 

w{T) = i^^^^^^^^^^f' =■■ A, w{r) > A, forO<r<T, (3.25) 

since 

w{T) = A> [e^^fTl^/", (3.26) 
because 1 < g < 5. So by ()3.H) and ()3.25|1 . we have 

f{r,w{r)) > 0, < r < T, 

and thus 

w'{r) < 0, for < r < T. (3.27) 
Remark 6. The existence of such a T follows here by a similar argument as in p!7] . 
The next lemma gives an upper bound for T, where the constant b is given by ()2.1(i|l . 

Lemma 3.2. Assume that 1 < q < 5 (6/5 < 7 < 2j and 

aS'{r)e''^z' > -^z, (3.28) 

for < r < n/^/b, where 

4r) = (3.29) 

then 

T < n/Vb. (3.30) 

Proof. It is easy to verify z{r) defined in ()3.29|) is the solution of the following initial 
value problem of second order linear equation 

z" + -z' + bz = 0, z(0) = 1, z'{0) = 0, (3.31) 
r 

By (HTT^ and (ICTj) . we have 

j^2gas{r)^^/^ - w'z)]' = r^[f{r, w)z - bze'^^w + a5'(r)e"^z'u;]. (3.32) 
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We prove ()3.30p by contradiction. Suppose T > ii/Vb, in view of ()3.1|) . ()3.26|) and 
f)3.27|) . we have 

w(r) 

> e--^(«;(T))'^-i - ^ - &e"^ = ^- (3-33) 
Integrating ()3.32j) over the interval [0, vr/ \/b], noticing the fact z{'k/ Vb) = 0, we obtain 

Jo y w z ) 

In view of (KTI^ . we have w(^) > w{T) > if T > ^. Moreover from (3.29) 
z'{^) < 0. Thus the left hand side of ()3.34|) is negative. On the other hand, by 
()3.28|) and ()3.33p . we can see that the right hand side of ()3.34p is positive. This is a 
contradiction, and thus ()3.30|) is proved. □ 

In order to prove Parts 2) and 3) in Theorem 2.1, we need a few lemmas. First, in 
view of ()3.1|1 and ()3.2|1 . there exists a positive constant B such that 

\e-^sir)p^^^ u;)! < 5, for < w < A, (3.35) 



e-"^^'^\g{r, w) - aS'{r)F{r, w)\<B, hi < w < A, (3.36) 

where A = w(T) is defined in ()3.25|1 . g{r, w) is defined in fl3.7j) . In the following, we 
denote 

Q :=Q{T) =w'iT), (3.37) 
where T = T{p) is defined in fl3.25|) . By ()3.27|) . we have 

g < 0. (3.38) 

Lemma 3.3. If the entropy S satisfies the hypothesis in Lemma 3.2, and 

45, A ^ VIA , iB AAB, 1 

and 

< w(r) < A = u;(T), for T < r < T + L, (3.40) 

for any L satisfying 

OA 

0<L< min{ — , 1}, (3.41) 
16 



then 

w'(r) < -i^, for T < r < T + L, (3.42) 
V2 

provided suPq<;^<;^_|_^ l'S"(r)| is sufficiently small. 

— Vh 

Proof. First, by ()3.3j) . we have 

+ e-"^F(r, w)y = _ 2 + «^'^'(^) ^/2 ^ g-a^Mr _ c,5'(r)F(r, u;)], (3.43) 

2 r 

where (^(r, w) is given by (jSIZI). Notice that T is bounded by n/Vb (cf. (2.30)), where 
b is given by ()2.16p . so 

[r,T + l] C [0,^ + 1]- (3.44) 
This, together with ()3.4ip . imphes 

[T, T + L]C [0,^ + 1]. (3.45) 

Therefore, if suPq<j,<^^^ l'S"(r)| is sufficiently small, then |rS"(r)| is also small for 
T <r <T + L. Thus 

< 2 + arS\r) < 3, (3.46) 
for T < r < T + L. It follows from (nOHll and (nOHll that, 

— + e-"^F(r, w) 

<9l + e-"^(^)F(T, u;(T)) + f [e-''^^'\g{t, w) - aS'{t)F{t, w)]dt, (3.47) 

for T < r < T + L. Using dSSSl), (!3:3fi|) and (j3l7jl . we obtain 

w'^ g2 g2 2BA 

_<^ + 2B + BL<^ + 2B + — (3.48) 

for T < r < T + L. On the other hand, we have from (3.36) and (3.43) that 

!f^+e-"^MF(r,w7) > ^+e-^(^)F(T,w7(r))- F l±^^^(w'itfdt-BL, (3.49) 
2 2 t 

for T < r < T + L. By (3.35), (3.46) and (3.49), we obtain 

\>\-^B-BL- \{w\tfdt 
> !^_25(1 + -)- / ^(u;'(t)2rft, (3.50) 



2 ' 2' T 
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for T < r < T + L. So, if L satisfies (3.41), by virtue of (3.48), we have 



2 >:|--2i?(l + ^)-y^ -^(uj'itrdt (61/ (3.41)) 



> ? - 2Bil + 4) - I + + 1# (3.48)) 



for T < r < T + L. Hence 

for T < r < T + L. Therefore, if (j3.39|l holds, then we have 

{w'{r)y > ^ (3.53) 

for T < r < T + L. This imphes w'{r) does not change sign for T < r < T + L. Since 
Q = w'{T) < (cf. (II2ZD), 

«;'(r) < -M (3.54) 
for T < r < T + L. □ 

Next, we have the following lemma. 
Lemma 3.4. If w{r) > for r G [T, T + L] where L satisfies (|J.^i| ), t/ien we /iave 



< w(r) < A, for T <r <T + L. (3.55) 

Proof. Since ti'(T) = A and w'{T) = Q < as we showed before, then w'(r) < 0, and 
thus w{r) < y4 for r > T, (r — T) small. We prove ()3.55p by contradiction. If ()3.55j] 
were false, then there exists r2 G (T, T + L] such that 

w{r) < y4 for r G (T, r2), ti'(r2) = A. (3.56) 
Since w{T) = w{r2) = A and r2 > T, by RoUe's Theorem, we have 

w'{t) = 0, (3.57) 
for some r G (T, r2). This contradicts ()3.54j) . The proof of the lemma is complete. □ 
The following lemma is a generalization of a result in |17] . 
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Lemma 3.5. Assume the entropy S satisfies the hypothesis in Lemma U~B . Let T{p) = 
T be the point defined in 25\j . If 

w'(T)T — oo, as p ^ +00, (3.58) 

then R{p, a) < +00 ifp is sufficiently large and suPq<^<j^_,_^ is sufficiently small. 

— — \/b 

Proof. First, in view of ()3.30|) . ()3.58p implies 

Q{T) = w'{T) —00, as p ^ +00. (3.59) 

Hence, there exists po > such that (3.59) holds for every p > pq. Now, for any fixed 
p G [pq, +00), we show that there exists r* G [T, T + ^] such that 

w{r*) < 0, (3.60) 

and this implies R{p, cr) < +00. We prove ()3.60|) by contradiction. Suppose 

2A 

w{r) > for r e [T, T + —]. (3.61) 

IWI 

Then by Lemma [3.41 we have 

2A 

< w{r) < A, for r G [T, T + —]. (3.62) 

IWI 

Thus, we can apply Lemma 13.31 with L = T + ^ to obtain 

y;'(r) < -M, {otT <r <T + ^. (3.63) 

' ' ~ V2 ~ ~ \Q\ 

Therefore, since w{T) = A, we have 



2A 



This contradicts 1)3.611) . The proof of the lemma is complete. □ 

For the case 1 < g < 3 (4/3 < 7 < 2), we have the following result. This is the same 
as Part 2) in Theorem 2.1. Proposition 1 Assume 1 < q < 3 (A/3 < 7 < 2^ and 
a > 0. If the entropy S{r) satisfies condition \3. 2^) in Lemma W^ then 

R{p, a) < +00, (3.65) 

for sufficiently large p, provzfiefi supo<^<i+^ 1*5" (r)| is sufficiently small. 

— \fb 
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Proof. To prove this proposition, it suffices to verify ()3.58|) in Lemma [3.51 This follows 
by the following argument. Set 

w'{T) = Q. 

By ()3.13j) . we have 

_^aS{T)j.2Q ^ f ^2j^^^ w{r))dr. (3.66) 







We estimate QT^ as follows. For < r < T, by ()3.27|) . we have 

w{r) > w{T). (3.67) 
Therefore, by ()3.1|) and ()3.25|) . we obtain 

/(r, w{r)) = e-°^(")w;«(r) - a 

> e-"^w''(T) - a 

> (i±lM^ > 0, (3.68) 

5 — g 

for < r < T, since 1 < g < 3. On the other hand, for ri defined in (3.14), since 
2p/3 > A = w{T) if p is sufficiently large, then (3.25) and ()3.27p imply 

T > n, (3.69) 

for p large. Hence, it follows from ()3.66|) - ()3.69|) that 

_^aS{T)r^2Q y f ' ^2j^^^ w{r))dr. (3.70) 
Jq 

By (3.14), we have 

w{r) > w{ri) > 2p/3, for < r < ri. (3-71) 

Thus, by (KT^ . 

f{r, w{r)) > e"^{3p/2y -a> canst ■ p", for < r < n. (3.72) 
if p is sufficiently large. Therefore, we obtain, by ()3.15|) . ()3.70|) and ()3.72|) . that 

_ga5(T)y2g > ^^^^^ . ^q^S ^> ^^^^^ . _ (^3 yg^. 

This implies T^Q — oo as p ^ +00 if g < 3. Condition ()3.58p is thus verified in view 
of in Lemma EH □ 
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For the case 1 < g < 5 (6/5 < 7 < 2), we have the following result. This is the same 
as Part 3) in Theorem 2.1. 



Proposition 2. Assume 1 < q < 5 (6/5 < 7 < 2^ and a > 0. If the entropy S{r) 
satisfies the hypothesis in Proposition 1 and 



IT 



S'{r) < 0, for < r < 



where b is given by \2. 1 6|) , then 



R{p, a) < +00, 



for sufficiently large p, provided supo<j.<i+^ is sufficiently small. 



(3.74) 



(3.75) 



Proof. First, in view of ()3.14|) and ()3.25|) . we have 

w(ri) > w(T), 

if p is sufficiently large. This, together with ()3.27p . implies 

ri < T. 

Once again, by ()3.25|) . we obtain 

w{r) > w{T), 
for ri < r < T. It follows from ()3.6|) that 

G(T) = G(ri) + [ 3r^ ( Pir, w) /(r, w)w H — rg(r, w)] dr. 

in V 6 3 J 

By dSH), dsn, dSIII), we have 

F(r, w) - ^/(r, w)w + lrg{r, w) 
D 3 



1 1 arS'(r) 



q+1 6 3(g + l) 



e-^'^wMy - -a}w(r) 
6 



for ri < r < T. If suPq<^<;i^^ l'S"(r)| is small, in view of ()3.30|) . we have 

1 1 -S'M^IM _l|>o. 



q + 1 6 3(g+ 1) 2 Vg + 1 6 
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(3.76) 



(3.77) 



(3.78) 



(3.79) 



(3.80) 



(3.81) 



for r < T and 1 < g < 5. Moreover, by (jT^ . (jSZZHl), and (jCTj) . we obtain 



1 1 
F(r, ly) - -/(r, w)w + -r5((r, w) 



e-'''{w{T)r--a\w{r)>(}, 



_q + l 6 

for ri < r < T. This, together with ()3.16p and ()3.79p . imphes 

G{T) > G{n) > C3p('-'')/^ 

By (2.5) and (jSSSl), we have 



(3.82) 



(3.83) 



2 



{[TQY + {TyF{T,w{T)) + 



^ [TQ] + ^ r^e'^'^^^S'{r){w'{r)fdr 

2 2 Jq 



(3.84) 



where Q = w'(T). So, if 5"(r) < for < r < we have, in view of (|3.3U|) . 

'-^[TQf + {TfFiTMT)) + —\ [TQ] 

i(5-q)/2 



Since T < ^ and 1 < g < 5, we have 

TQ — >• — oo, as p — »• +00. 
By virtue of Lemma 13. 5^ Proposition 2 is proved. 



(3.85) 

(3.86) 
□ 



Now we prove Part 4) of Theorem 2.1. For the solutions w(r, a) =: w{r) in 2) 

1/(7-1) 



and 3) of Theorem 2.1, let p 



7-1 



W 



M{r) 



; p is the density function (cf. (2.3)). Set 

47rr^p(r)(ir, 



the mass in the ball i?,.(0). We calculate each term in (2.5) as follows, by virtue 
of the fact that w and S are spherically symmetric and f2 is a constant. First, for 
r = Y/xf~Ta?|~Tx|, we have 

Vw = Wr{xi/r, x^jr^ x^/r), VS" = S'r(xi/r, X2/'r, Xs/r), (3.87) 



and thus 



ae WrSr + e (Wrr + -^r) 

r 

e"Vr + e"^(- + a5'(r))w,. 
r 



(3. 
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By (2.3) and (2.6), we have 

= AirGp. (3.89) 
Substituting ()3.88j) and ()3.89j) into (2.5), and noticing Q is a constant, we obtain 

e^^Wrr + e"^(- + aS'{r))wr + AuGp - 2^'^ = 0. (3.90) 



r 



Thus 



(r^e^'^w'ir))' + r^AnGp - 2^^) = 0, (3.91) 
where G is the Newtonian gravitational constant. Therefore 

r2e°V(r) = — GM{r), (3.92) 

for r < R{p, a). By the proofs of Part 2) and 3) of Theorem 2.1, we have w'{r) < 
for r < i?(p, a). Thus M(r) > for r < R{p, a). The estimate fTI^ follows 

immediately from ()2.22|) . 



4 Proofs of Theorems |2H and |2l3 

In this section, we consider the isentropic case S = constant. Theorem 2.1 can be 
applied in this case because the conditions imposed on the entropy S in Theorem 2.1 
are automatically satisfied when S is constant. Without loss of generality, throughout 
this section, we assume 

S = (4.1) 
for convenience. We prove Theorems 12.21 and 12.31 as follows. First, we define 

f(w) = - a, w>0, (4.2) 

and 

F(w) = aw, w>0 (4.3) 

q + I 

We still use w{r,p, a) to denote the solution of the problem 

w"(r) + -w'(r) + fiw) = 0, (4.4) 
r 

and 

w'{0) = 0,w{0) =p> 0. (4.5) 

Let R{p, a) be defined as in (2.13); R{p, a) is the first point at which w is 0. Then 
we have the following proposition. 
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Proposition 3. If a > and 1 < q <5, (6/5 < j < 2), then 

R(p^ a) = +00, if <p<p, (4.6) 

where 

,.(5fll±i))./,.(i52^p-.. (4.7) 

- — q 07 — 



Proof. It is easy to verify that i^^^^^Y^'^ is the unique positive zero of F{w) — 
If R{p, a) < +00, we let R{p, a) = R > 0. Then G{R) > 0, where the function G is 
defined as in (2.5) with S = 0. On the other hand < w{r) < p for < r < R. If 
p<p = (5£li±9l)i/g^ then {F{w) - ^^)(r) < for < r < i?. Hence, (El implies 
G{R) < 0. This is a contradiction. □ 

Proof of Part 1) of Theorem 2.2 
For fixed cr > 0, let 

Pq = mi{p : R{p, a) < +00}. (4.8) 
By Part 3) of Theorem 2.1 and Proposition 3, we know that 

P<Po< +00, 

where p is given by ()4.7j) . If we can show 

R{po, a) < cx), (4.9) 

then Part 1) in Theorem 2.2 will be proved. Now (|4.9j) can be shown by the following 
argument. Write ()4.4j) as a first order system 

2 

w' = V, V = — V + cr — w'^, (4-10) 
r 

with initial condition 

w{^) =po>0, viO) = 0. (4.11) 
Define the Hamiltonian H{w, v) by 

^2 g+l 

H{w,v) = - + ^^-aw. (4.12) 

Then H' = —^v^, so that H decreases on orbits of ()4.10|) . This implies that the solution 
of ()4.1U|) - ()4.1ip must have v bounded from below in the region S = {w > 0, f < 0}. 
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This solution also cannot exit S via w = at some v < for some R > 0; otherwise, 
since v'{R) = —j^v{R) + cr > 0, the solution crosses the line w = transversally, so by 
continuity, there would be a neighborhood of po on the w-axis, such that for p E N, 
the orbit of ()4.1()|1 satisfying w{0) = p, v{0) = 0, would also exit S* at a point near 
w = 0, v{R) < 0. This would contradict the definition of po- Similarly the po orbit 
cannot exit S via some point {w, 0), with < w < Pq. Thus the po orbit exits 5* via 
w = 0, V = 0, so (2.31) holds, and this proves Part 1) of Theorem 1.2. 

Remark 7. It is not hard to show that every solution of ()4.1()j) satisfying w{0) = p, 
f (0) = 0, tends to the rest point {w = a^^^ , i; = 0) as r ^ oo. 



Proof of Part 2) of Theorem 

First, we show that the radius of the star increases with the increasing angular 
velocity. 

We use Wj(r) (i = 1, 2) to denote the solution to the following initial value problem 

^"(r) + f<(r) + w'i-a, = 0, r> 0, 

w,iO)=p, wm = o- 

We assume 

ai > 0-2 > 0, (4.14) 

and we want to show that 

Rip, (Ti) > Rip, (T2). (4.15) 

(Note that each of these are finite, by what we have already shown in Part 2) of the 
theorem) . 

From ()4.13|) . we have 

wlir) = \ f s^ai- wlis)) ds, ^ = 1, 2, (4.16) 
Jo 

It is easy to verify, using L'Hospital's rule, that 

lua-w[ir) = lia,~p'^). (4.17) 
It follows from and (HTfll that 

lim Wi"ir) = -ip" - ai)/3. (4.18) 

r— >0+ 
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Since ai > (T2, ()4.18|) implies 

w"i{0+) > w"2{0+). (4.19) 

This, together with the fact Wi{0) = ^2(0) = p and w[{0) = ^2(0) — 0; leads to 

Wi{r) > W2{r), for small r > 0. (4.20) 

We shall show ()4.15p by contradiction. If ()4.15|) were false, then there exists an tq 
< ro < R{p, (Ti) < R{p, (72) such that 

wi{r) > 1^2 (r), for < r < tq, and wi{ro) = 1^2 (ro). (4-21) 

Let 

y(r) = wi{r) — W2{r). 

We then have from (4.13), that 

y"(r) + V(r) + yB{r) - {a^ - a^) = 0. (4.22) 

Here ^ 

B{r) = q I {\wi + (1 - X)w2Y~^ {r)d\. (4.23) 
Jo 

From (|i?T^ and dOT)), we have 

1/(0) = ?/'(0) = 0, y{r) > for, < r < Tq and ?/(ro) = 0. (4.24) 

Multiplying ()4.22|) by y' and integrating the resulting equation over the interval [0, ro], 
we get, since |/'(0) = 0, 

^ Jo ^ 

rro pro 

+ / B{r)yy'{r)dr -I ((Xi - (J2)y'{r)dr = 0. (4.25) 
Jo Jo 

With the help of ()4.24j) and integration by parts, we get 

fro rro ^2 i rro 

B{r)yy'{r)dr = B{r){—)'dr = — / B\r)y^{r)dr (4.26) 



Jo 2 2 Jo 



and 



pro 

/ (ai - a2)i/'(r)dr = 0. (4.27) 
Jo 
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substituting and KTi\i into (jOSl), we have 

yyyo)) ^ / f(^/(^))2^^_^^ / ^-B'{r))y\r)dr = 0. (4.28) 



2 Jo ^ 2 

By the definition of B{r) (see fl4.23|) ). we have 

B' {r) = q{q - I) f (Au;i + (1 - A)u;2)''~^(AK(r) + (1 - A)w^(r))dA. (4.29) 
Jo 

Since w-(r) < for < r < R{p, ai) (^ = 1, 2) (P), we thus have 

B'{r) < 0, for < r < Tq, (4.30) 

when 7 < 2, i.e., g > 1. Hence , each term in ()4.28|) must be zero. This contradicts 
dOl), and proves (gUS)). 

Now we show that the radius of the star decreases with the increasing central density, 
i.e., 

i?(pi,fx) > i?(p2,o-), (4.31) 

if Po < Pi < P2, and a > 0. For this purpose, let w{r,p,a) be the solution of the 
following initial value problem 



(4.32) 



w"{r) + f w'(r) + 10" -a = 0, 
w{0) =p, w'{0) = 0, 
for p > po- We use the following rescaling, 

A = rp'^''-^^/\ e{X) = w/p ; (4.33) 

then 6{X) is the solution of the following initial value problem: 

^AA + f ^^A + - ^ = 0, A > 0, 

^(0) = 1, ^'(0) = 0. 

The first zero of 6 depends only on the parameter ^. We use A(^) to denote this first 
zero. Then by ()4.33|) . we have 

i?(p,a)=p(i-«)/2A(^), (4.34) 

for p > Po- Similar to the argument in the proof of ()4.15|) . we can show that A(^) 
increases with the parameter ^. Thus, for the fixed a, it decreases with p for p > po- 
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This, together with ()4.34|) and the fact q > I, imphes ()4.31|) . This completes the proof 
of Part 2) of Theorem 1.2. 



Remark 8. The above scahng argument also works for non-rotating star, i.e., the case 
when (7 = 0. For the non-rotating star, the radius of the star, R{p, 0), is always finite 
for p>0, if6/5<7<2 (see For the non-rotating star, ()4.34|1 becomes 

i?(p,0) = Aop(^-^)/^ q=^, (4.35) 

7-1 

where Aq is the first zero of the function 6{X), which is the solution of the following 
initial value problem, 

^AA + f + = 0, A > 0, 

^(0) = 1, ^'(0) = 0. 

From (|4.35p we see that the radius of a non- rotating star is proportional to ^(^^9)/^. 

We now turn to the proof of Theorem 2.3. In (j4.34j) . since A(^) increases with the 
parameter ^, as we have already shown, we have 



0" . , , cr 



< M-^), (4.36) 
P Po 

for p > Pq. We apply ()4.34|) to the case of p = po and obtain, 

RiPo,cr)=pt'^^'X{^). (4.37) 
Po 

Now R{po, cr) < +00, from Part 2) of Theorem 2.1, and po is a positive constant deter- 
mined completely by 7 and a. Thus A(^) is a positive constant determined also only 
by 7 and a. We set this positive constant as C in ()4.3fij) . then ()2.31|) follows. Since 
7 < 2, ()2.3H1 implies ()2.32j) . The proof of Theorem 2.2 is complete. 

The following lemma will be useful in proving Theorem 12.41 

Lemma 4.1. For any cr > 0, there exists Rq > such that \2.2!^] with the data 
w'(0) = and w{R) = has a positive solution in the ball Br{0) for all R < Rq, if 
l<q<5, f6/5 < 7 < 2). 

Proof. Fix cr > 0, and write R{p, o") = R{p). By Theorem 12.21 Part (1), there exists 
Po > such that R{p) < +00 for any p > pq. Let Rq = R{po)- For any < R < Rq, 
it follows from ()2.32|) that, there exits pi > po such that R{pi) < R. Now R{po) > 
R > R{pi), and R{p) is a continuous function of p (actually, it is differentiable; cf. 
[T7]). Thus, by the intermediate value theorem, there exists p G {po,Pi) such that 
R{p) =R. □ 
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5 Proofs of Theorems 12. 4L 2.5 and 2.6 

In this section, we consider the case for the general domains and variable angular 
velocity. In Theorems 2.4 and 2.5, we assume the entropy function S{x) is constant. 
We set S* = for convenience. Substitute this in ()2.5p . we get 

Aw + w'^ - 2n{ri){n{ri) + r]n\r])) = 0, (5.1) 

We look for the solution of ()5.1|1 satisfying the following condition 

w{x) >0,x e D, w{x) = 0, X edD. (5.2) 



In order to prove Theorem 12. 4| we need a comparison lemma, which can be found in 
[Zj or • Before we state this lemma, we give the following definitions of weak sub- 
and supersolutions of the problem 1)5.11) and ()5.2j) . For notational convenience, define 
the function f{x, w) by 

/(x, w) = w'^- 2n{ri{x))[n{ri{x)) + r/(x)fi'(r/(x)]. 



Definition, (a) w E H^{D) is called a weak supersolution of problem h5. 1\) and h5. g)) if 



Vw-Vvdx> I f{x,w)vdx (5.3) 
Jd 

for each v G Cq{D), v > a.e., where the subscript zero denotes v = on dD. 
Similarly 

(b) WE H^{D) is called a weak subsolution of problem \5.1]) and i5.^) if 

/ Ww_-Vvdx< / f{x,w)vdx (5.4) 
Jd Jd 

for each v G Cq{D), v >0 

The following lemma is well-known; cf [7j or |16j . 

Lemma 5.1. Assume there exists a weak supersolution w and a weak subsolution w_ of 
( 15. ip and \5.^) satisfying 

U7 < 0, w >Q on dD in the trace sense, and w<w a.e. in D. (5.5) 

Then there exists a weak solution w G Hq{D) of l[5.1]} and i5.^} . such that 

w < w < w a.e. in D. (5.6) 
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Remark 9. By the smoothness assumptions of dD and ^{f]), the usual regularity argu- 
ments (see [lOj) show that a weak positive solution of ()5.1|) and ()5.2|) must be a classical 
solution. 

Let D be a bounded domain with smooth boundary, and assume ^(r/) G C^^^{D) fl 
C{D), for some 5, < 5 < 1. Define /5 by 

(3 =:max\2n{r]){n{r]) + r]n'{r])\, (5.7) 

where r] = rj{x) is as defined in Section 1. Let i?/jj(0) be the biggest ball contained in 
D. Consider the positive solution of the equation 

u"{r) + -u\r) + u« - /? = 0, (5.8) 
r 

with data 

n'(0) = 0, u{r) > for < r < u{Ri) = 0. (5.9) 

By Lemma [4. 11 if6/5<7<2, then there exists a positive constant R2 depending only 
on P and q such that (j5.8|) and (j5.9j) has a solution if < Ri < R2. We require 

D C BrM- (5.10) 

This implies that (j5.8|) and (|5.9p has a positive solution. We denote this solution by 
u{x) = u{\x\) {x G i?R^(0)). A subsolution to ()5.1|) and ()5.2|) can now be constructed. 
First define w{x) by 

f u{\x\) for X G 5/j,(0) 
w x) = < ^ ^ (5.11) 

^ I for X G D-Br^{0). ^ ^ 

Lemma 5.2. // Condition A in Section 1 holds, then w{x) defined by i5.ll]) is a weak 
subsolution to i5.1]) and i5.^) . 

Proof. First, by the definition of w , we have 
(Vw ■ Vf — f{x, w)v) dx 

D 

(Vu ■ Vf — /(x, u)v)dx 



+ / 2Q{ri{x))[Q{ri{x)) +r]{x)Q'{r]{x)]v{x)dx (5.12) 

for each v G Cq{D), v > 0. By Condition A in Section 1, we have 

n{ri{x))[n{ri{x)) + ri{x)n'{ri{x)] < 0, (5.13) 
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for X E D — i?/j^(0). On the other hand, since u satisfies fj5.8|) in the ball Br^^O), we 
have by the divergence theorem, 

Vm ■ Vvdx 

Auvdx + / —vdS 

r r du 

< / f{x,u)vdx+ / -T^vdS, (5.14) 

where u is the unit outer normal vector. Since u'{R) < , we have |^ < on dBji^{0), 
and since u{x) = w{x) for x G i?j?j(0), ()5.13|1 implies 



\/u-\/vdx< / f{x,w)-vdx, (5.15) 
so that 

{Vw-Vv-f{x,w)v)dx<0, (5.16) 

Jd 

and this proves Lemma 15.21 . □ 

We next construct a supersolution to ()5.1|) and ()5.2|) . Let -R2 be the positive number 
as in ()5.10|1 such that D C i?ii;2(0) and p be any positive number such that 

p > maxw;(a;), (5-17) 

x&D 

where w is the subsolution which we have already constructed. We consider the fol- 
lowing boundary value problem 

Aw + w'^ + P = 0, for a; G ^^^(O), (5.18) 

and 

w\dBii^{o) = P- (5.19) 

By a result in j^, there exists a unique positive solution w{x) to problem ()5.18|1 and 
()5.19|1 . Furthermore, by the maximum principle , we have 

w{x) > p, for X G Br^{0). (5.20) 

Therefore, by the choice of p (see f|5.17p . we have 

w{x) > w(x), for X E D. (5-21) 
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It is easy to verify that w satisfies ()5.3p . i.e., w{x) is a supersolution of ()5.1|) and ()5.2|) . 
Then fl5.5|) follows from ()5.2H) and the properties of w_ and w. Hence, by Lemma 4.1, 
problem (|5.1|) and (j5.2j) has a solution w satisfying w_{x) < w{x) < w{x) for x (z D. 
Inequality ()2.34j) can be proved in the same way as ()2.28|1 . Now we prove w{x) > for 
X E D. Since w{x) > w{x) for x E D, and w{x) > if x E i?R^(0), it suffices to show 
w{x) > a.s X E D — B^-^{0). This can be proved by the strong maximum principle 
(|10j). using Condition A in Section 1. This completes the proof of Theorem 12.41 

We now prove Theorem 12.51 For this, we recall that the solution w{x) of ()5.H1 and 
(|5.2p is constructed by the following iteration (cf. [7j or [E])- Set Wo{x) = w{x), and 
then define Wk {k = 0, 1, 2, ■ ■ ■ ) inductively to be the unique weak solution of the 
linear boundary-value problem 

-Awk+i + Cwk+i = f{wk) + Cwk, in D 
Wk+i = 0, on dD, 

where C = max^g^ \f'{w)\{x). Then it can be shown (cf. (THj) 

w = Wq < Wi < ■ ■ ■ < Wk < ■ ■ ■ < w , a.e. in D. 

Set w{x) = limk^oo uik{x); then w{x) is the desired solution. Since w and w are axi- 
symmetric, if D is an axi-symmetric domain, it is not hard to show each Wk {k=0, 1, 
2, ....) is axi-symmetric; i.e. depends only on rj and z. Indeed, if 7] = \/xl + x\, and 
z = X3, then the Laplacian transforms to 

T] 

So wi satisfies the equation 

-{dnnWi + ^dr^wi + d^zWi) + Cwi = f{w) + Cw, (5.22) 

together with the boundary condition wi = on dD. Since the right hand side of ()5.22p 
depends only on rj and z, and since this Dirichlet problem, being linear has a unique 
solution, it follows that wi is axi-symmetric. Similarly, each Wk is axi-symmetric. Thus 
w{x) is axi-symmetric. This completes the proof of Theorem 12.51 

We finally prove Theorem 12. 6| the proof makes use of the celebrated Pohozaev iden- 
tity (see HH). 
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Proof of Theorem \2.(A 

Multiply (2.5) by x ■ Vw and w respectively, and integrate the resulting equations 
over D. After some manipulation (details can be found in [14]), one obtains the fol- 
lowing " Pohozaev-like " identity 



' D 



4 



q + l' 



aKe 



-aS 



D 



q + l 



-w 



''+^]{x-VS)dx 



+ / (riA'ir]) + 



5A{7]) 



D 



wdx+ / -e''^\Vw\'^{x-u)dS = 0, 

JdD 2 



(5.23) 



where A{rj) = 2Q{ri){Q{ri) + rjVL'ij]), and v is the unit outer normal vector. If 7 < 6/5 
(resp. 7 < 6/5), then g > 5 (resp. g > 5) and thus \ — > (resp. | — > 0) . 
Thus, if r]A'{r]) + ^ > (resp. r]A'{ri) + ^ > 0), x ■ z/ > and a; ■ V5 > 0, (HT:^ 
implies w{x) = for x E D. This completes the proof of Theorem 12.61 

6 Further Discussion. 

In this section, we give further discussions on the difference between our formulation 
of the problem and that in [2] and [11]. As mentioned in Section 1, for the isentropic 
case, i. e., 5* = constant, system (1.3) is considered in [2] and [11] in all of i?^ space, 
with the prescribed total mass M. In [2] and [11], the function $ is given by 

$(x) = -G [ T^^dy. (6.1) 

V \x - y\ 

This choice is natural when one considers the following equation A$ = AttGp in all of 
R^, because any bounded solution $ for this equation in differs from (j6.H) only by 
a constant. 

In our formulation, we consider system (1.3) in a bounded domain D, the domain 
of a star. The boundary condition in our formulation is p = 0, to match the exterior 
vacuum continuously. In this formulation, the function $(a;) could be different from 
the following formula 

^x) = -G [ j^dy, (6.2) 

Jd f - y\ 
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due to the boundary effects. We illustrate this by considering the isentropic case, i. e., 
S — constant (without loss of generality, S is assumed to be zero). In this case, the 
first equation in (1.3) reduces to 

V(w + * - J) = 0, (6.3) 

where w is given by (2.3) (with S — 0), J is given by (1.4). From ( 6.3), we have 

w + ^-J^C, (6.4) 

where C is a constant. On the boundary dD, since 7 > 1, the boundary condition 
p = implies w = on the boundary dD. This, together with (6.4), implies that 

$(x) = C + J, for X e dD. (6.5) 

The solution of the problem A$ = AnGp with the boundary condition (6.5) could be 
different form formula (6.2). The difference is a harmonic function which counts the 

boundary effects. 

When the domain D is ball and the entropy S* is a constant, we give the existence 
results and study some physical properties of the solutions. Historically, incompress- 
ible fluid model of rotating stars are studied extensively. For the incompressible fluid 
model, density is a positive constant inside a star and zero outside the star, so there 
is a discontinuity of density across the boundary of a star. The problem for the in- 
compressible fluid model of Newtonian rotating stars is to determine the boundaries 
of stars with the prescribed angular velocity (or angular momentum) and total mass. 
For such problems, explicit solutions are found (e. g. Maclaurin spheroids and Jacobi 
ellipsoids) . Those spheroids and ellipsoids are close to balls if the angular velocities are 
small. For the compressible fluid model, the problem is quite different since the distri- 
bution of density must be determined. For the compressible fluid model of a rotating 
star, if the domain D is not a ball, it seems very difficult to obtain some interesting 
properties of solutions like what we obtain in Theorems 2.1-2.3 for the case when D is 
a ball. When D shghtly differs from a ball, for example, an ellipsoid with eccentricity 
close to one, it is reasonable to expect the solutions obtained in this paper for the case 
when D is a ball give good approximations. 

Some non-existence results are proved in [11] for large angular velocity Q for the 
problem formulated in [2]. This means, in order to ensure existence of solutions for 
the problem formulated in [2]; the angular velocity cannot exceed a critical value. This 
critical value is given in term of total mass M, which is prescribed in the formulation 
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of the problem in [2] and [11]. Our existence theorem, Theorem 2.1, is vahd for any 
angular velocity Q. This difference is due to the different formulations of problems 
mentioned above. In our formulation, we do not prescribe the total mass, instead we 
prescribe the domain. The total mass thus depends on the angular velocity and the 
domain. We also give a lower bound of the total mass (cf. (2.22)) in terms of angular 
velocity. 
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